x = (0:20);
y = nbinpdf(x,3,0.5);
plot(x,y,'+')
